Seismic Electric Signals and 1/f "noise" in natural time 
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By making use of the concept of natural time, a simple model is proposed which exhibits the 1//" 
behavior with a close to unity. The properties of the model are compared to those of the Seismic 
Electric Signals (SES) activities that have been found to obey the ubiquitous 1//" behavior with 
a « 1. This comparison, which is made by using the most recent SES data (that were followed 
by three magnitude 6.0-class earthquakes), reveals certain similarities, but the following important 
difference is found: The model suggests that the entropy S- under time reversal becomes larger 
compared to the entropy S in forward time, thus disagreeing with the experimental SES results 
which show that S may be either smaller or larger than S- . This might be due to the fact that SES 
activities exhibit critical dynamics, while the model cannot capture all the characteristics of such 
dynamics. 

PACS numbers: 05.40.-a, 05.45.Tp, 91.30.Dk, 89.75.-k 
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I. INTRODUCTION 

Among the different features that characterize com- 
plex physical systems, the most ubiquitous is the pres- 
ence of 1//" noise in fluctuating physical variablesfl]. 
This means that the Fourier power spectrum S{f) of 
fluctuations scales with frequency / as S{f) ~ I//"- 
The power-law behavior often persists over several or- 
ders of magnitude with cutoffs present at both high and 
low frequencies. Typical values of the exponent a ap- 
proximately range between 0.8 and 4 (e.g., see Ref.[3| 
and references therein), but in a loose terminology all 
these systems are said to exhibit 1/f "noise". Such a 
"noise" is found in a large variety of systems, e.g., con- 
densed matter systerns (for example, an excellent review 
can be found in Ref.^l), freeway traffic [1, d, 01 , granular 
flowQ, DNA sequence Q , heartbeat @ , ionic current fluc- 
tuations in membrane channels [l ^, r iver discharge HtII, 
the number of stocks traded daily chaotic quantum 
systemslllflllil, [il], the light of quasars p7j. human 
cognition] 18 1 and coordination . burst errors in com- 
munication systems [2(|, electrical measurements (2l| . the 
electric noise in carbon nanotubesjl^l and in nanoparticle 
fllmspsj. the occurrence of earthquakes (2^ etc. In some 
of these systems, the exponent a was reported to be very 
close to 1, but good quality data supporting such a value 
exist in a few of them[3]. As a flrst example, we refer 
to the voltage fluctuations when current flows through a 
resistor ^2^. As a second example we mention the case 
of Seismic Electric Signals (SES) activities which are 
transient low frequericy (< IHz) electric signals observed 
before earthquakes Hi, 113, [H, [H, H^, IM [H, [H, [H, 
since they are emitted when the stress in the focal re- 
gion reaches a critical value before the failure [ssl. [36| . 
These electric signals, for strong earthquakes with mag- 
nitude 6.5 or larger, are also accompanied by detectable 



TABLE I: The values of 5", 
presented in Fig[S] 



Ki, S- for the electric signals 



Date recorded 


S 




S- 


Feb 8, 2007 
Apr 23, 2007 
Apr 24, 2007 

Nov 7, 2007 


0.067±0.007 
0.07f±0.005 
0.072±0.003 
0.070±0.005 


0.074±0.007 
0.069±0.003 
0.067±0.003 
0.065±0.005 


0.079±0.007 
0.066±0.005 
0.069±0.003 
0.070±0.005 
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magnetic field variations [371, l38|, l39|, |40|. Actually, the 
analysis of the original time series of the SES activities 
have been shown to obey a 1/ f-behavior[4li[4^. 

The 1//° behavior has been well understood on the 
basis of dynamic scaling observed at equilibrium crit- 
ical points where the power-law correlations in time 
stem from the infinite-range correlations in space (see 
Ref. and references therein). Most of the observations 
mentioned above, however, refer to nonequilibrium phe- 
nomena for which -despite some challenging theoretical 
attempts [46l. l47l l48l . l49l |- possible generic mechanisms 
leading to scale invariant fluctuations have not yet been 
identified. In other words, despite its ubiquity, there is 
no yet universal explanation about the phenomenon of 
the 1//° behavior. Opinions have been expressed (e.g., 
see Ref ■ [Tsj ) that it does not arise as a consequence of 
particular physical interactions, but it is a generic mani- 
festation of complex systems. 

It has been recently shownUl M, HH, M, HI, [H, [H, 
m, 113) HI) US US im that novel dynamic features hidden 
behind the time series of complex systems can emerge 
if we analyze them in terms of a newly introduced time 
domain, termed natural time x (see below). It seems that 
this analysis enables the study of the dynamic evolution 
of a complex system and identifies when the system enters 
a critical stage. Natural time domain is optimal [6^ for 
enhancing the signal's localization in the time frequency 
space, which conforms to the desire to reduce uncertainty 
and extract signal information as much as possible. In a 
time series comprising N events, the natural time Xk = 



2 




FIG. 1: (color online) (a):Example of the evolution of e„ versus the number of renewals n, i.e., in natural time. An exponential 
PDF has been considered for the selection of rjn (see the text), (b): The Fourier power spectrum of (a); the (green) solid line 
corresponds to 1/f and was drawn as a guide to the eye. (c): The DFA of (a) that exhibits an exponent an fa very close to unitj, 
as expected from (b). (d):Properties of the distribution of e„. The average value (e„) (plus) and the variance {(e„ — (en)) ) 
(crosses) as a function of n. The straight solid line depicts ln(n) and was drawn for the sake of reader's convenience. 



k/N serves as an index [4 ll ISOl l5l| for the occurrence 
of the fc-th event. The evolution of the pair (yfc, Qfc) is 

studiedlSi, SH, M, M, M, M, M,M,mM,^, 

where Qk denotes a quantity proportional to the energy 
released in the fc-th event. 



The scope of the present paper is twofold. First, a sim- 
ple model is proposed (Section |TT| which, in the frame of 
natural time, leads to 1//° behavior with an exponent 
a close to unity. Second, the properties of this model in 
natural time are compared to those of the SES activities 
in Section [TTTl This comparison is carried out by making 
use of the most recent experimental data of SES activ- 
ities observed in Greece during the last several months. 
Section [IV] presents the conclusions. Two Appendices 
are also provided, the first of which refers to the earth- 
quakes that followed the SES activities presented here. 
Appendix B clarifies that an early model proposed by 
Voss[44] differs essentially from the one presented here. 



II. THE MODEL PROPOSED 



A. Description of the model 



Here, we present a simple competitive evolution model 
which results, when analyzed in natural time, to 1//° 
"noise" with a very close to unity. Let us consider the car- 
dinality e„ of the family of sets En of successive extrema 
obtained from a given probability distribution function 
(PDF); Eo equals to the empty set. Each En is obtained 
by following the procedure described below for n times. 
Select a random number ri„ from a given PDF (here, we 
use the exponential PDF, i.e., p{rin) = exp(— 77„)) and 
compare it with all the members of E,i_i. In order to 
construct the set E,i, we disregard from the set En-i 
all its members that are smaller than rjn and further- 
more include rjn. Thus, En 7^ for all n> and En is 
a finite set of real numbers whose members are always 
larger or equal to 77„. Moreover, min[En] > min[E„_i] 
and max[En] > max[En_i]. The cardinality = |En| 
of these sets, which may be considered as equivalent to 
the dimensionality of the thresholds distribution in the 
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FIG. 2; (color online) Results from 10* runs of the model 
presented in Fig[T] (a) the average power spectrum, (b) De- 
trended Fluctuation Analyses of order / (DFA-Z)[43]. The 
black solid line in (a) corresponds to 1// spectrum and was 
drawn as a guide to the eye. For the same reason in (b), the 
black solid lines correspond to udfa = 1. In (b), the colored 
solid lines correspond to the least square fit of the average 
Fbfa-i, depicted by symbols of the same color; the numbers 
in parentheses denote the standard deviation of aoFA-i ob- 
tained from the 10* runs of the model. The various Fofa-i 
have been displaced vertically for the sake of clarity. 



coherent noise model (e.g. see Ref. 62] and references 
therein), if considered as time-series with respect to the 
natoraZ number n (see Fig|IJa)) exhibits 1//° noise with 
a very close to unity, see FiglTJb) . This very simple model 
whose evolution is depicted in Fig[lja), leads to a De- 
trended Fluctuation Analysis 9] (DFA) exponent an pa 
close to unity, see FigHJc), being compatible with the 
1/f power spectrum depicted in Fig[T]^b). The mathe- 
matical model described above corresponds to an asymp- 
totically non-stationary process, since (e„) cx Inn with a 
variance ((e„ — (e„)) ) (x Inn (see FiglTJd)). Thus, in 
simple words, the present model suggests that the cardi- 
nality e„ of the family of sets En of successive extrema 
exhibits 1//° behavior when considered as time-series 
with respect to the natural (time) number n. We note 
that a connection between 1//° noise and extreme value 
statistics has been established and proposed as providing 
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FIG. 3: The r;„ values arranged in sites (bins) according to 
their value increasing from left to right. (a)The six(=3!) 
equally probable outcomes after the selection of 3 random 
numbers by the same PDF. Actually, the sample space is (in 
one to one correspondence to) the permutations of 3 objects, 
(b) The 24(=4!) equally probable outcomes after the selection 
of 4 random numbers by the same PDF. Again, the sample 
space is (in one to one correspondence to) the permutations of 
4 objects. For the reader's convenience, in each outcome, the 
corresponding en-value (n = 3 or 4) is written. An inspection 
of (b), shows that p(e4 = 1) = l/4,p(e4 = 2) = 11/24, p(e4 = 
3) = 1/4 and p(e4 = 4) = 1/24. 



a new angle at the generic aspect of the phenomena [4a| . 
Furthermore, in the frame of a formal similarity between 
the discrete spectrum of quantum systems and a dis- 
crete time series [3| the following striking similarity is 
noticed: The fact that a « 1 together with the behavior 
((fn ~ (^n))^) cx Inn of the present model is reminiscent 
of the power law exponent and the statistic in chaotic 
quantum systems [T3, [IB| . 

In order to check the stability of the results of Fig[l] we 
present in Figl^Ja) the average power spectrum obtained 
from 10* runs of the model. A sharp 1/f behavior is 
observed. Moreover, in Fig[2fb), we present the results 
of the corresponding average values of FjjpA-i obtained 
from DFA of various orders I (i;^ when detrending with a 
polynomial of order I, see Ref . pSj ) . Figure [2tb) indicates 
that aoFA-i is close to unity. 



B. Analytical properties 

We now discuss an analytical procedure which clarifies 
some properties of the model. In order to find analyti- 
cally the distribution of the probabilities p(en), one has 
simply to consider the possible outcomes when drawing 
n random numbers rjn. Since the selection is made by 
a means of a PDF, all these numbers are different from 
each other, thus -when sorted they- are equivalent to n 
points (sites) lying on the real axis. The value of e„ varies 
as {r]n} permutate along these n sites independently from 
the PDF used in the calculation. Thus, a detailed study 
of the permutation group of n objects can lead to an ex- 
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FIG. 4: (color online) (a)The probabilities p(en — m) as a 
function of m for various n. The bar chart corresponds to the 
exact p(e3 = m) whereas the continuous lines to the Cornish- 
Fisher approximation of Ea. (|15|l . The latter approximation 
converges very rapidly to the true p(en = rn), see for example 
n = 3. This fact enables the calculation of p{en ~ m) for 
very large n, for which the recursive relation of Eq.([lJ would 
accumulate significant round-off errors, (b) The probabilities 
p(i'2'"-i ~ j) a-s a function of j for various m (see Appendix 
B). They are clearly skewnessless, i.e., symmetric around their 
mean. Here, as in Ref . [4^ . 6-sided dices (fc — 6) were consid- 
ered. For the reader's convenience p{v2m-i — j) for m — 3 
and m = 12 have been drawn with a different style. 
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FIG. 5: Four electric signals recorded at PAT (sampling rate 
fexp^l sample/sec) on February 8, 2007(a), April 23, 2007(b), 
April 24, 2QQ7(c) and November 7, 2007(d). 



act solution of the model. It is well known, however, that 
the number of the elements of this group is n! and this 
explains why we preferred to use the numerical calcula- 
tion shown in FigU] Some exact results obtained by this 
method are the following: (ei) = 1; (£2) = 1 + 1/2, since 
Pie2 = 1) = p{e2 = 2) ^ 1/2; (eg) = 1 + 1/2 + 1/3, since 
Pies - 1) = l/3,p(e3 = 2) = 1/2 and p(e3 = 3) = 1/6; 
(£4) = 1 + 1/2+1/3 + 1/4 (see Fig^. Figure [3] analyzes 
the results for n = 3 (FiglSlJa)) and n = 4 (Figlltb)). 
One can see that the probability p{en = ™) equals to the 
sum of the n possible outcomes as rjn moves from the left 
to right in the n columns of FiglJl In each column, the 



probability to have at the end e„ = m is just equal to 
the probability to keep m — 1 numbers from the numbers 
already drawn that are larger than 77,1 . This results in 



p(e„ = m) = - p{ek = m - 1) (1) 

n ^ — ' 

k—m~l 

(cf. p(eo = 0) = 1). 

Equation ([1]) enables us to calculate the characteristic 
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FIG. 6: (color online) (a): The electric signal depicted in 
FigEl^c) (April 24,2007) in normalized units (i.e., by subtract- 
ing the mean value and dividing the results by the standard 
deviation) along with its dichotomous representation which is 
marked by the dotted (blue) line, (b): How the signal in (a) 
is read in natural time. 



function (see p.928 of Ref.flil) 

n 

/„(A) = (exp (Ae„)) - ^ e^Xe„ = m). (2) 



FIG. 7: (color online) Evolution of the parameters of ki, S 
and S- as a function of n, when e„ are analyzed in the nat- 
ural time domain. The thick lines correspond to the average 
value of Ki, S and S-, found by 10* runs of the model. The 
thinner lines correspond to the ±one standard deviation con- 
fidence intervals. For the reader's convenience, the green and 
black horizontal lines show the values Ku and Su of ki and S, 
respectively, that correspond to a "uniform" distribution. 



Since /i(A) — e'^. we find that Eq.® -upon considering 
Eq.6.1.22 of Ref.jill- results in 



/n(A) 



1 r(e^ 



r(e^) 



(6) 



where T{x) is the gamma function. Now, the mean and 
all the central moments fii = ((e„ — (e„))') of the distri- 
bution of p{en = TTi) can be obtained by virtue of the 
cumulant theorem (see p.928 of Ref . [6^ ) : 



7n—l k—m-'l 



Subtracting now Eq.Q from Eq.Q, we finally get 



/„+l(A) ^ ^^^fnW. 
71+1 



(4) 



(5) 



Indeed, by substituting p{en — m) in Jn{^)^ we obtain 




{en) = 




n n—1 

nU\)=Y.^^^ E P{^k=m-\), (3) 

m— 1 k—m — 1 


At2 = ((en 


-(Of> = 


^ln/,.(A) 


whereas by substituting p(en+i — Tn) in /n+i(A), we find 


1^3 = ((fn 






n n — 1 

(n + l)/„+i(A) - ^ e^'" J2 p{ek=m-l) + e^h{n). 




pi4 - = 





A=0 



A=0 



A=0 



A=0 



(7) 
(8) 
(9) 
(10) 



Substituting Eq.® into Eqs.iO to Uni) and using the 
properties of the polygamma functions (i.e., the n-th or- 
der logarithmic derivatives of the gamma function, see 
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FIG. 8: (color online) Evolution of the parameters of ki, S 
and S- as a function of n, when Vn (see Appendix B) are an- 
alyzed in the natural time domain. Since the model produces 
1/f noise time series with lengths n — 2™, the calculation 
was performed at such n (i.e., n = 10™, m = 4, 5, . . . 16) in- 
dicated by the points marked. The thick lines connect the 
corresponding points of the average value of /ti, S and S-, 
found by 10* runs of the model. The thinner lines connect 
the points corresponding to the ±one standard deviation con- 
fidence intervals. For the reader's convenience, the green and 
black horizontal lines show the values and Su of ki and 
S, respectively, that correspond to a "uniform" distribution. 
Here, as in Ref.'43|, 6-sided dices (k — 6) were considered. 



p.260 of Ref.lH), we obtain 
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(11) 

(12) 
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Equations (fTTj) to ((T4|) enable us to calculate the mean, 
standard deviation cr(= y/J^), skewness 71 — [i^jcr' and 
kurtosis 72 = /i4/o'^ — 3 as a function of n. Using now 
the Cornish-Fisher (CF) expansion treated in Ref. !69j , we 
obtain the following continuous approximation to p(e„ = 
m) 



X exp ■ 



PCF(£71 

~2 



/27r 



2l,~„ + 2i(i2e;,2 
3 36 ^ 



7) 



y (^n' - 1) - - 36;0 + ||(4e;3 



where — (cn — (£«))/"■• Equation (fT5|) . although be- 
ing a continuous approximation to the point probabili- 
ties p(e„ = m), rapidly converges to the latter, see for 
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FIG. 9: The map shows the areas A,B. The star indicates the 
epicenter of the strong 6.0 EQ that occurred on March 25, 
2007 in Kefallonia. 

example the comparison of the exact p(e3 = m) and the 
corresponding pcf(£3)/c in FigHlE^)- An inspection of 
this figure, which depicts the probabilities p(e„ = m) up 
to n = lO"^, reveals that even for large n, the probability 
p(e„ — m) remains non-Gaussian (cf. even a.t n — 10^, 
we obtain from Eqs.(IT2]) to ^ 71 = 0.2154 ^ with 
72 = 0.0459 > 0). 



III. COMPARISON OF THE MODEL WITH 
THE SES PHYSICAL PROPERTIES IN 
NATURAL TIME 

For dichotomous signals, which is frequently the case 
of SES activities, the quantity Qk mentioned in Section I 
stands for the duration of the fc-th pulse. The normalized 
power spectrum Il{uj) = |<i>([jj)p was introduced ^4l]. [sol. 
ISlj . where 



N 



^(uj) = Pk exp 



fe=i 



N 



(16) 



and Pk = Qkl Qn, ^ = 27r(/); 4> stands for the natu- 

ral frequency. The continuous function <^{uj) should not 
be confused with the usual discrete Fourier transform, 
which considers only its values at (/i = 0, 1, 2, . . .. In nat- 
e~j^ral_):iiRe analysis [Sa. l4ll. [sol. Islj . the properties of n((jj) 
or n((p) are studied for natural frequencies (p less than 0.5, 
?i| this range of 0, n(a;) or 11(0) reduces to a charac- 
f ti(Yfi,^tion for the probability distribution pk in the 
contex-t of probability theory. When the system enters 
the critical stage, the following relation holds [ill. Isol lE^ : 

18 6cosw 12 sin w 
5^ 
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FIG. 10: (color online) (a) The normalized power spectrum(red) 11(0) of the seismicity as it evolves event by event (whose 
date and time (UT) of occurrence are written in each panel) after the initiation of the SES activity on February 8, 2007. The 
excerpt presented here refers to the period 19 to 25 March, 2007 and corresponds to the area B, Mthres = 3.2. In each case 
only the spectrum in the area G [0,0.5] is depicted (separated by the vertical dotted lines), whereas the 11(0) of Eq.(2) is 
depicted by blue color. The minor horizontal ticks for 4> are marked every 0.1. (b), (c) Evolution of the parameters (D), ki, 
S and S- after the initiation of the SES activity on February 8, 2007 for the areas B {Mthres = 3.2) and A{Mthres = 3.2), 
respectively, until just before the 6.0 EQ. 



For uj^Q, Eq.inD leads to[3a |4l|, l50| 

n(tj) w 1 - o.07cj2 

which reflects [131 ^^^^ the variance of x is given by 



(x'> - ixf - 0.07, 



where (/(x)) = T,k=iPkfiXk)- The entropy S in the 
natural time-domain is defined as[50l l53l| 

5" = (xlnx) - (x> ln(x>, 

which depends on the sequential order of events 5^, UHl . 
It exhibits [5^ concavity, positivity, Leschejcl, [67| stabil- 
ity, and for SES activities (critical dynamics) its value 
is smallerjH HI] than the value Su{= In 2/2 - 1/4 w 

0. 0966)of a "uniform" (u) distribution (as defined in 
Refs. [53, m, m, m, HHI, e.g. when all pk are equal or 
Qk are positive independent and identically distributed 
random variables of finite variance. In this case, ki and S 
are designated Ku{— 1/12) and Su, respectively.). Thus, 
S < Su- The same holds for the value of the entropy 
obtained [H, [5§] upon considering the time reversal T, 

1. e., Tpk = pN-k+i, which is labelled by S-. 

In summary, the SES activities, when analyzed in nat- 
ural time exhibit infinitely ranged temporal correlations 



and obey the conditions |59|, |60|: 
Kl = 0.07 

and 

S,S- < Su. 



(18) 



(19) 



We first present in subsection llll Al the most recent ex- 
perimental results on Seismic Electric Signals activities 
and their properties are compared, in subsection IIII Bl 
with those of the model proposed as well as with a dis- 
crete model proposed j44| by R.F.Voss(see Appendix B). 
The classification, in advance, of these SES activities, on 
the basis of relations ^8]) and p9|) has been verified by 
the occurrence of three magnitude 6.0-class earthquakes 
(see Appendix A). 



A. The recent electric field data 

Figure [5] depicts the original time series of four elec- 
trical disturbances that have been recently recorded on: 
(a) February 8, 2007, (b) April 23, 2007, (c) April 24, 
2007 and (d) November 7, 2007 at a measuring sta- 
tion termed Patras (PAT) located at ~ 160km west of 
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FIG. 11; The map shows the areas C and D. The star shows 
the epicenter of the strong 5.8 EQ at 18:09:11 on June 29, 
2007. 

Athens. All these four recent electric signals were ana- 
lyzed in natural time. For example, if we read in natural 
time the signal on April 24, 2007 (FigEfc)) -the dichoto- 
mous representation of which is marked by the dotted 
(blue) line in FiglB^a)- we find the natural time represen- 
tation of FiglUb) the analysis of which leads to the values 
Ki = 0.067±0.003, S = 0.072±0.003, S- = 0.069±0.003. 
The relevant results of all the four signals arc compiled 
in Tableland found to be consistent with the conditions 
(fT8|) and (fT9|) , thus they can be classified as SES activities 
(for their subsequent seismicity as well as for more recent 
SES activities see Appendix A). An inspection of Table |T] 
shows that the S value is more or less comparable to that 
of 5*-, but experimental uncertainty does not allow any 
conclusion which of them is larger. Note that in several 
former examples[56], the data analysis also showed that 
the S value may cither be smaller or larger than S-. 



B. Comparison of the SES properties with those of 
the model proposed 

We now turn to investigate whether the parameters 
Ki, S and S- deduced from the 1/f model of Section HIl 
are consistent to those resulted from the analysis of the 
SES activities observed. Figure [7] summarizes the results 
of 10"* runs of the model which, for moderate sizes of 
n, seems to obey more or less the conditions ([TBI) and 
(fT9|) . In particular, for n < 10^ (which is frequently 
the number of pulses of the SES activities observed in 
field experiments), Fig[7| shows that ki is close to 0.070, 
S < Su and (in most cases) S- < Su- A closer inspection 
of FiglTl however, reveals the following incompatibility of 
the model with the experimental results: For n < 10^, 



the model clearly suggests that S- > S, thus disagreeing 
with the experimental data which show, as mentioned 
above, that S may either be smaller or larger than 5_. 
The origin of this incompatibility has not yet been fully 
understood. It might be due to the fact that SES ac- 
tivities exhibit critical dynamics, while the model cannot 
capture all the characteristics of such dynamics. 

Since the model proposed here might be considered as 
reminiscent of an early 1/f model proposed by R.F.Voss 
(see Ref. [iS] ) , which is also a discrete one, we also present 
in FiglH] the corresponding results of 10"* runs of that 
model. (The details of the model are given in Appendix 
B). A comparison of Figs [7] and [5] for n < 10^ reveals that 
the results of the two models differ essentially. For ex- 
ample, in the Voss model, S- is larger than Su while the 
opposite holds for the model proposed here. Moreover, 
when comparing the results of the Voss model with those 
of the SES activities we also find considerable differences. 
For example, the ni value deduced from the Voss model 
is (on the average) larger than k„, thus differing consid- 
erably from the value ki « 0.070 of the SES activities. 
This, as mentioned above, is comparable to the one de- 
duced from the model proposed here and explains why we 
focused on this paper on the comparison of that model 
-among the variety of models suggested to date for the 
explanation of the 1/f behavior- with the experimental 
values obtained from the analysis of the SES activities 



IV. CONCLUSIONS 

In summary, using the newly introduced concept of 
natural time: (a) A simple model is proposed that ex- 
hibits 1//" behavior with a close to unity, (b) Electric 
signals, recorded during the last few months in Greece, 
are classified as SES activities since they exhibit infinitely 
ranged temporal correlations. Actually, three magnitude 
6.0 class earthquakes already occurred in Greece (see Ap- 
pendix A), (c) For sizes n comparable to those of the 
SES activities measured in the field experiments (i.e., 
n < 10^), the model proposed here leads to values of 
the parameters ki (« 0.070) and S,S^ (< Su) that are 
consistent with those deduced from the SES activities 
analysis. Despite of this fact, however, the model results 
in S- values that are almost always larger than those of 
5, while the observed SES activities result in S values 
that may be either larger or smaller than S-. This dis- 
crepancy might be due to the inability of the model to 
capture the characteristics of critical dynamics which is 
exhibited by SES activities. 



APPENDIX A: WHAT HAPPENED AFTER THE 
SES ACTIVITIES DEPICTED IN FIGJl 

We clarify that, during the last decade, preseismic 
information[68] based on SES activities is issued only 
when the magnitude of the strongest EQ of the im- 
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FIG. 12: (color online) (a),(b) and (c),(d) depict the evolution of the parameters (D), ki, S and S- after the initiation of the 
SES activity on April 24, 2007 for the areas C and D, respectively (for two magnitude thresholds in each area), until 03:40:15 
UT on June 25, 2007. 
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FIG. 13: (color onhne) The map shows the areas N^EioI 
and A'^Iy I iSlo o in which the seismicity was studied [43 after 
the SES activity recorded at PAT on 7 November, 2007. The 
solid dots stand for the sites at which electric field variations 
are continuously monitored with a sampling frequency IHz. 



pending EQ activity is estimated -by means of the SES 
amplitude[20, [11, [H, [sO]- to be comparable to 6.0 
units or larger [3g. Here, in the first two subsections, we 
explain what happened after the SES activity at PAT 
on February 8, 2007 (see FiglHJa)) and on April 23 and 
24, 2007 (see Figsl5ljb),(c)). The relevant analysis of seis- 
micity after the SES activity on November 7, 2007, which 
was still in progress during the initial submission of 
this paper on 23 November, 2007 and hence completed 
afterwards, is presented in Subsection 3. 



1. What happened after the SES activity of 
February 8, 2007 

According to the National Observatory of Athens, 
NOA (the seismic data of which will be used here), a 
strong earthquake (EQ) with magnitude 6.0-units oc- 
curred at Kefallonia area, i.e., 38.34°N 20.42°E, at 13:57 
UT on March 25, 2007. We show below that the occur- 
rence time of this strong EQ can be estimated by follow- 
ing the procedure described in Refs. [36l. [sol [stI. IsqI. [60j . 

We study how the seismicity evolved after the record- 
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FIG. 14: (color online) The evolution of the parameters (D), 
Ki, S and S- for the period: since the initial submission of 
this paper until just before the 6.6 earthquake on 6 January, 
2008. (a), (b) correspond to the two areas depicted in Fig |13l 
while (c) to the dark shaded area in Fig[T5j Mthrea ~ 3.1. 



ing of the SES activity at PAT on February 8, 2007, by 
considering either the area AiiVgy g SIq q or its smaller 
area BiA^l^ li^lo.o (see FigH]). If we set the natural time 
for seismicity zero at the initiation of the SES activity on 
February 8, 2007, we form time series of seismic events 
in natural time for various time windows as the num- 




ia"E 19"E 20"E 21 "E 22"E 23"E 24"E 25"E 26"E 27"E 2a"E 29"E 30"E 

FIG. 15: (color online) The map shows the epicenters of the 
5.3 earthquake on 30 December, 2007 (small star) and the 
6.6 earthquake on 6 January, 2008 (large star). The dark 
shaded area depicts the part of the larger area A'^ff;! iJlo.'o 
adjacent to the epicenter of the 6.6 earthquake, whose the 
four parameters (D), ki, S and S- fuUfilled the conditions for 
a true coincidence just one day before the major earthquake 
occurrence (see Fig[l4jc)). 



ber N of consecutive (small) EQs increases. We then 
compute the normalized power spectrum in natural time 
n(0) for each of the time windows. Excerpt of these 
results, which refers to the values deduced during the pe- 
riod from 20:53:59 UT on March 19, 2007, to 11:56:30 UT 
on 25 March, 2007, is depicted in red in FigHOja). This 
figure corresponds to the area B with magnitude thresh- 
old (hereafter referring to the local magnitude ML or the 
'duration' magnitude MD) Mthres = 3.2. In the same 
figure, we plot in blue the power spectrum obeying the 
relation (2) which holds, as mentioned, when the system 
enters the critical stage. The date and the time of the 
occurrence of each small earthquake (with magnitude ex- 
ceeding (or equal to) the aforementioned threshold) that 
occurred in area B, is also written in red in each panel. 
An inspection of this figure reveals that the red crosses 
approach the blue line as TV increases and a coincidence 
occurs at the last small event which had a magnitude 
3.2 and occurred at 11:56:30 UT on March 25, 2007, i.e., 
just two hours before the strong 6.0 EQ. To ensure that 
this coincidence is a true one (see also below) we also cal- 
culate the evolution of the quantities ki,S and S- and 
the results are depicted in Fig. [TUfb) and [TUTc) for the 
same magnitude thresholds for each of the areas B and 
A, respectively. 

The conditions for a coincidence to be considered as 
true are the following (e.g., see Refs.[36l[50l.[57l.[59ll60t): 
First, the 'average' distance {D) between the empirical 
and the theoretical 11(0) (i.e., the red crosses and the blue 
line, respectively, in Fig fTOT a)) should be smaller than 
10~^. See Fig. [Wb).(c) where we plot {D) versus the 
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FIG. 17: The electric signal recorded on 9 October, 2008 at 
PAT in a fashion similar to that of Fig[5] 
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FIG. 16: (color online) The electric signals recorded on 10 
January, 2008 at PAT (a) and on 14 January, 2008 at PIR 
(b), in a fashion similar to that of Fig|6ja) 



2. What happened after the SES activities of April 
23 and 24, 2007 

We investigate the seismicity after the aforementioned 
SES activities depicted in FigslSjb) andOJc). The in- 
vestigation is made in the areas C: ^37 g £^20 2 ^^"^ 
-^38.1 ^20.2 fsee Fig fTTj) . Starting the computation of seis- 
micity from the initiation of the SES activity on April 24, 
2007 (which, between the two SES activities depicted in 
FigsEI^b) and[5]^c), has the higher actual amphtude), we 
obtain the resuhs depicted in Figs[T2fa),(b) and flW c^.fd) 
for the areas C and D, respectively, for Mthres = 3.0 and 
Mthres = 3.2. An inspection of the parameters (D), Ki,S 
and S- reveals that they exhibited a true coincidence (as 
discussed above) around June 25, 2007, i.e., around four 
days before the 5.8 EQ that occurred at 18:09:11 UT on 
June 29, 2007, with an epicenter at 39.3°N 20.3°E (shown 
by a star in Fig fTTj) . 



conventional time for the aforementioned two areas B 
and A, respectively. Second, in the examples observed to 
date, a few events before the coincidence leading to the 
strong EQ, the evolving 11(0) has been found to approach 
that of the relation (2), i.e., the blue one in Fig fTUf a) , 
from below (cf. this reflects that during this approach the 
Ki-value decreases as the number of events increases). In 
addition, both values S and S- should be smaller than 
Su at the coincidence. Finally, since the process con- 
cerned is self-similar ( critical dynamics) , the time of the 
occurrence of the (true) coincidence should not change, 
in principle, upon changing the surrounding area (and 
the magnitude threshold Mthres) used in the calculation. 
Note that in Fig. [TOTb). upon the occurrence of the afore- 
mentioned last small event at 11:56:30 UT of March 25, 
2007, in area B the (D) value becomes smaller than 10^^. 
The same was found to hold for the area A, see FigfTCTc). 



3. Study of the seismicity after the SES activity on 
November 7, 2007 

This study, which as mentioned above was still in 
progress upon the initial submission of this paper, was 
made by investigating the seismicity in the area B of 
FigEl i.e., ^37'|£^2Qn ^^^^ ^ the larger area, i.e., 
^37.1 ^20.0 (see Ref.;45]), which are shown in FigUl The 
parameters ki, S, and (D) computed during the sub- 
sequent period for Mthres = 3.1 are depicted in FigsfT4ya) 
and [T^b) for the smaller and larger area, respectively. 
An inspection of these two figures reveals that the con- 
ditions for a true coincidence (see subsection A.l) were 
obeyed upon the occurrence of the first event early in the 
morning (i.e., at 03:25 UT) on 30 December, 2007 with an 
epicenter at 37.8°N 20.2°E and magnitude 3.9. Almost 
three hours later, i.e., at 06:42 UT, a strong earthquake 
of magnitude 5.3 occurred at 37.6°N 20.9°E marked with 
the small star in Fig |15l In addition, and quite interest- 
ingly, the four parameters ki, 5, S- and (D) during the 
next few days continued to fulfill the conditions for a true 
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coincidence, as it is evident from a closer inspection of 
Figsimja) andlTUb). Actually, at 05:14 UT on 6 Jan- 
uary, 2008, a major magnitude 6.6 earthquake occurred, 
which was felt not only all over Greece but also in adja- 
cent countries, e.g., southern Italy and western Turkey. 
Its epicenter, marked with a large star in FigllSl was lo- 
cated at 37.1°N 22.8°E, i.e., only around 50 km to the 
south of the larger area studied since the initial submis- 
sion of this paper. Interestingly, in the part of the latter 
area adjacent to the epicenter (which is shaded in FigfTS]) 
the four parameters ki, S, and (D) reached the con- 
ditions for a true coincidence just in the morning (i.e., 
at 04:24 UT) of 5 January, 2008 upon the occurrence of 
a magnitude 3.6 earthquake at 38.4°N 22.0°E (it corre- 
sponds to the last point in FigfT^c) where (D) becomes 
smaller than 1%). 

Finally, we note that after the initial submission of this 
paper, two additional SES activities have been recorded 
as follows (see FigHH): One SES activity at PAT on 10 
January, 2008 and another one on 14 January, 2008 at 
the station FIR located in western Greece, see FigfT51(cf. 
The configuration of the measuring dipoles in the latter 
station is described in detail in the EPAPS document 
of Ref.(6C|). Their subsequent scismicities are currently 
studied along the lines explained above considering the 
evolving seismicity in the following areas: Concerning 
the former SES activity at PAT the areas depicted in 
Figini while for the one at PIR on 14 January, 2008, the 
subsequent seismicity is studied in the area B of FiglH] 
as well as in the larger area -/VH qE'Iq q and in the one 
surr ounding the epicenter^ (36°N 23°E). 

We now offer some comments on the classification of 
the aforementioned electric signals of Fig [11] as SES activ- 
ities. Concerning the signal on 14 January, 2008, which 
is of clear dichotomous nature, the analysis is made by 
considering that Qk stands for the duration of fc-th pulse, 
as mentioned in Section I, and the following parameters 
are obtained: ki = 0.070 ± 0.005, S = 0.086 ± 0.003, 
S- = 0.070 ± 0.005, which obey the conditions UHl) 
and (fT9|) for the classification of the signal as SES ac- 
tivity. Furthermore, note the S- is smaller than that 
of S, which is not compatible with the model proposed 
in Section II, thus strengthening the point mentioned in 
Section III.B as well as in the Conclusions (Section IV) 
that the model does not seem to capture the character- 
istics of critical dynamics exhibited by SES activities. 
We now turn to the signal recorded at PAT on 10 Jan- 
uary, 2008, the feature of which is not clearly dichoto- 
mous since it consists of three main pulses that seem 
to overlap (Note that, in general, if pulses of very short 
duration exist, the calculation of Qk should necessarily 
consider the characteristics of the low pass filters used in 
our measurements [36|, i4Q!i- This has been considered in 
drawing the dichotomous representation -marked by the 
dotted (blue) line- in FigfrBFa')). Its analysis leads to the 
following values ki = 0.070 ± 0.010, S = 0.050 ± 0.010, 
S- = 0.060 ± 0.010. These values, which are different 
from those deduced from the analysis of the SES activ- 



ity on 7 November, 2007, also obey the conditions (fT8|) 
and (fT9|) . At this point, we clarify that the optimality of 
natural time domain for enhancing the signal's localiza- 
tion in the time- frequency space was shown ^6^ without 
assuming that Qk stands for the pulse duration, but it 
was noted that in general Qk is a quantity proportional 
to the corresponding energy released in the fc-th event 
(estimated by means of the time integration of the signal 
exceeding the level of the irrelevant background noise). 



4. Additional data 

Meanwhile, the earthquakes that followed the SES ac- 
tivities: (i) depicted in Fig[T6] and (ii) recorded on 29 
February - 2 March, 2008 lat er at P IR (see N.V.Sarlis, 
E.S. Skordas, M.S. Lazaridou, [arXi7:0802.3329v4), have 
been described in N.V.Sarlis, E.S. Skordas, M.S. Lazari- 
dou, Proc. Jpn. Acad. B 84, 331, 2008. In addition, 
a strong electrical disturbance has been recorded at 9 
October, 2008 fFig fTT]) . having an amplitude comparable 
to that in FigO^d). The relevant analysis of the sub- 
sequent seismicity in natural time is carried out in the 
area N^^fEfg^. Such an analysis, when considering the 
NOA seismic data published until December 6, 2008, we 
find the results depicted in Figs. UM a) and UMh) for 
magnitude thresholds 3.0 and 3.2, respectively. An in- 
spection of these figures shows that Prob( ki) exhibits a 
local maximum at ki w 0.070, thus obeying EqlTHl The 
spatial invariance of this result is currently investigated 
in order to check whether the critical point has been ac- 
tually approached. 



APPENDIX B: COMPARISON OF THE 
PRESENT MODEL WITH A MODEL 
SUGGESTED BY R.F.VOSS 

A model producing 1/f noise, suggested by Richard 
F. Voss, was presented in Ref . p3| . This model as- 
sumes some number m of fc-sided dices, i.e., the outcome 
of a rolling the l-th dice is equally distributed among 
the values di = 1,2,... A:. To each dice I, one relates 
the l-th digit of the binary representation of a number 
n — 0,1, ... , 2™ — 1. The procedure to generate the 1/f 
noise starts (n = 0) by rolling all dices and assign their 
sum to vo{= '^i)- second step (n = 1), one 

reroUs only the dice associated with the least significant 
digit of the binary representation, and the new roll d'l 
is summed with the previous rolls of all other dices so 
that vi = d'l + 'Ylii=2 ^i- T^^^ procedure continues up to 
n — 2"^ ~ 1, each time rolling only the dices associated 
with the digits that change when considering the binary 
representation of n compared to n—1. This model results 
in time-series w„ of length 2™ whose spectrum is close to 
1 /f , and their values Vn are obviously distributed among 
the integers m and mk. The time series v„ clearly differ 
from e„ apart from the fact that they both have integer 
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FIG. 18: (a) Prob(Ki) versus Ki of the seismicity, for Mthres = 3.0, subsequent to the SES activity at PAT on October 9, 2008. 
The last event corresponds to the Ml = 3.0 event that occurred at 00:43 UT on December 6, 2008 at 37.81°iV21.15°i;. (b)The 
same as in (a), but for Mthres = 3.2 
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values. To visualize this difference, we depict in Fig|4l^b) 
the distribution of p{v2m-i) for the case of 6-sided dices 
(fc = 6). It can be easily found, since upon consider- 
ing the binary representation of n = 2™~^ compared to 
n — 1 — 2™~^ — 1, all digits change and one reroUs all 
dices. Clearly, rolling all dices results in a symmetric 
(skewnessless) distribution for V2m-i (see Fig|4l[b)). In- 
deed, by considering the distribution of the sum of rolling 
m independent fc-sided dices and using the characteris- 
tic function method discussed in subsection lllBl one can 
find the following cumulants: 



Clearly, Eqs. (|BT1) to ((B4| for the distribution of Wa^-i 
differ from Eqs.pTj) to (fTi]) for the distribution of e„. 
Among their differences, the following two are the most 
striking: First, V2m-i is skewnessless (see Ea. (|B3p ) 
whereas that of e„ is not, and second the two distribu- 
tions have different signs in their kurtoses. 
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